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Abstract. As a model of a neutron halo nucleus we consider a neutron bound to an inert core by a zero 
range force. We study the breakup of this simple nucleus in the Coulomb field of a target nucleus. In the 
post-form DWBA (or, in our simple model CWBA ("Coulomb Wave Born Approximation")) an analytic 
solution for the T-matrix is known. We study limiting cases of this T-matrix. As it should be, we recover 
the Born approximation for weak Coulomb fields (i.e., for the relevant Coulomb parameters much smaller 
than 1). For strong Coulomb fields, high beam energies, and scattering to the forward region we find a 
result which is very similar to the Born result. It is only modified by a relative phase (close to 0) between 
the two terms and a prefactor (close to 1). A similar situation exists for bremsstrahlung emission. This 
formula can be related to the first order semiclassical treatment of the electromagnetic dissociation. Since 
our CWBA model contains the electromagnetic interaction between the core and the target nucleus to all 
orders, this means that higher order effects (including postacceleration effects) are small in the case of high 
beam energies and forward scattering. Our model also predicts a scaling behavior of the differential cross 
section, that is, different systems (with different binding energies, beam energies and scattering angles) 
show the same dependence on two variables x and y. 



PACS. 25.70.De Coulomb excitation 
and momentum distributions 



25.60.-t Reactions induced by unstable nuclei - 25.60.Gc Breakup 



1 Introduction 

Breakup processes in nucleus-nucleus collisions are com- 
plicated, in whatever way they are treated. They con- 
stitute at least a three-body problem, which is further 
complicated due to the long range Coulomb force. Exact 
treatments (like the Faddeev-approach) are therefore pro- 
hibitively cumbersome. On the other hand, many approx- 
imate schemes have been developed in the field of direct 
nuclear reactions, and these approaches have been used 
with considerable success Jj}. In this context we wish to 
investigate a realistic model for the Coulomb breakup of a 
neutron halo nucleus. With the operation of exotic beam 
facilities all over the world, these reactions (previously re- 
stricted essentially to deuteron induced reactions) have 
come into focus again. The Coulomb breakup of these 
nuclei is of interest also for nuclear astrophysics, since 
the breakup cross section can be related to the photo- 
dissociation cross section and to radiative capture reac- 
tions relevant for nuclear astrophysics g. 

An important benefit of the present model is that it 
can be solved analytically in the DWBA (CWBA) ap- 
proximation. Thus it constitutes an ideal "theoretical lab- 
oratory" to investigate the physics of breakup reactions, 
certain limiting cases and its relation to other models like 
the semiclassical approximation, which is mainly used in 



the interpretation of experiments. Especially the effect of 
postacceleration (to be explained in more detail below) 
can be studied in a unique way in this approach. 



2 Description of Theoretical Model 

We consider the breakup of a particle a = [c + ri) 
(deuteron, neutron-halo nucleus) consisting of a loosely 
bound neutral particle n and the core c (with charge Z c ) 
in the Coulomb field of a target nucleus with charge Z. 



c + n + Z. 



(1) 



As a further simplification the a = (c + n) system is as- 
sumed to be bound by a zero range force. The bound-state 
wave function of the system is given by 



K exp(— nr) 



(2) 



where the quantity k is related to the binding energy Sbind 
of the system by 
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In the post-form CWBA the T-matrix for the reaction 
Eq. (|l|) can be written as || 



T 



Vn, 



D J ^^(RJe-^-^^R), 



(4) 
(5) 




The 



) range constant'' 
initial state 

,(+) 



D given by D = 
is given by the incoming 

Coulomb wave function x<u' with momentum q a and the 
halo wave function <f>Q. The final state is given by the in- 
dependent motion of the core described by the outgoing 

Coulomb wave function Xq c ' in the Coulomb field of the 
target nucleus Z with asymptotic momentum q c and the 
free neutron with momentum q„, described by a plane 
wave. In these wave functions the Coulomb interaction is 
taken into account correctly to all orders. In our model 
there is no resonance structure in the c + n continuum. 
This is clearly a good assumption for the deuteron and 
also for other neutron halo systems. 

There exists another form of the T-matrix element, 
which is not equivalent to Eq. (||). It is called the "prior- 
form" Q. The final state is described by a cm. motion of 
the (c + n) system (as a Coulomb wave function) and a 
relative wave function of the unbound (c + n) system. 

The present "post-form" description, Eqs. (jj), ^ in- 
cludes the effects of "postacceleration" . "Postaccelera- 
tion" arises in a purely classical picture of the breakup pro- 
cess. This is nicely discussed in || (We show their figure 5 
here as our Fig. Fy). The nucleus a = (c + n) moves up the 
Coulomb potential, loosing the appropriate amount of ki- 
netic energy. At the "breakup point" (marked as "breakup 
occurs here", see Fig ||), this kinetic energy (minus the 
binding energy) is supposed to be shared among the frag- 
ments according to their mass ratio (assuming that the ve- 
locities of c and n are equal) . Running down the Coulomb 
barrier, the charged particle c alone (and not the neutron) 
gains back the Coulomb energy, resulting in its "postac- 
celeration" . Of course this picture is based on the purely 
classical interpretation of this process, and will be modi- 
fied in a quantal treatment, where such a "breakup point" 
does not exist. The correct semiclassical limit of the the- 
ory in this case can be found, e.g., in ||. A purely classi- 
cal formula for this postacceleration, where the "breakup 
point" corresponds to the distance of closest approach — 
i.e., b — r in Fig. [j] — is given in Postacceleration is 
clearly observed in low energy deuteron breakup, in the 
(fully quantal) theoretical calculations as well as in the 
corresponding experiments, see Fig. ||and also, e.g., @,|[- 

The formula Eq. (^) is also useful for the description of 
the Coulomb dissociation of halo nuclei at high beam ener- 
gies, see 1 10 . Within this theory postacceleration effects 
become negligibly small in the high energy region. This 
is seen in the numerical calculations pj and in the ana- 
lytical investigations to be described below. It can, e.g., 
be applied to 11 Be and 19 C Coulomb dissociation experi- 
ments [|ll],[l2| (We disregard here the importance of finite 
range effects). 



a=(c+n) 




photon 
absorbed 
here 



breakup 
occurs 

here 



Fig. 1. A schematic view of Coulomb-breakup, as adapted 
from Q. The distance from the target nucleus to the breakup 
point is denoted by r. 




Fig. 2. Comparison of calculations and measurement for the 
deuteron breakup coincidence cross section on 197 Au at Ed = 
12 MeV (Fig. 4 of fl). The postacceleration effect can clearly 
be seen, as the maximum of the proton energy (~ 7.5 MeV) is 
larger than the one of the neutron (~ 2.5 MeV). The experi- 
mental data are taken from fl9|. 



On the other hand the I s * order semiclassical Coulomb 
excitation theory was widely applied in the past years to 
the Coulomb dissociation of high energy neutron halo nu- 
clei, see, e.g., Q. The theory corresponds to the "prior 
form" , mentioned above. The question of hi ghe r order 
electromagnetic effects was studied recently in WM within 
this framework. These effects were found to be small, for 
zero range as well as finite range wave functions of the 
a = (c + n)-system. It seems interesting to note that 
postacceleration effects arises through higher order elec- 
tromagnetic effects in straight line semiclassical theories, 
see |l5j. Through the interference of I s * and 2 nd order am- 
plitudes even a "post-deceleration" can arise, as was seen 
in that paper. 

In this work we want to establish the relation between 
the apparently very different post-form CWBA and semi- 
classical theory. It was recently noticed [jl3| that in the 
limit of Coulomb parameters r\ a = Z c Ze 2 /hv a <C 1 (i.e. in 
the Born approximation) , where v a denotes the velocity of 
particle a (v a = £2a), both theories give the same result. 
Expanding the Coulomb wave functions up to first order 
in the Coulomb fields one finds 
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qint=q n +^x^ 



qcoui 
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Fig. 3. The two bremsstrahlung type of graphs, which de- 
scribe the Coulomb breakup in the Born approximation. Three- 
momentum conservation at each vertex determines the inter- 
mediate momenta qi nt . 
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Here / cou i = 2r] a q a / (qcoui) is the usual Coulomb ampli- 
tude with the "Coulomb push" q cou i = q a — (q c + qn), for 
further details see j}). The two terms in theparenthesis 
correspond to the two graphs shown in Fig. |. For small 
values of q cou i the two terms almost cancel and the ex- 
pansion in q cou i was found to be in agreement with the 
semiclassical result, see and below. 

We now show that this agreement is also true in the 
case of arbitrary values of rj a and rj c . The beam energy 
must be high (compared to the binding energy T'&md) and 
the two fragment need to be scattered into forward angles. 

This is reminiscent of the result in the theory of 
bremsstrahlung. Replacing the neutron by a photon the 
diagrams of Fig. || are identical to the bremsstrahlung in 
lowest order. In this case it was already noticed that the 
Born result remains valid for arbitrary values of r\ for high 
energies and small scattering angles [16 . Here we want to 
show that the same applies in this case. 

The T-matrix can be evaluated analytically in this 
model due to well known Nordsieck formula fl^l , see 
Eqs. (11)-(13) of [||. Using this formula one obtains the 
T-matrix Eq. (fj) in terms of a hypergeometric function F 
as well as its derivative F'. The argument of the hyperge- 
ometric function F (and F') is given by P,Iuj: 



C(A) 



(7) 



2 9coul(9a9c 



q a q c ) - 4(q cou iq Q + A(7 a )(q cou iq c - Xq c ) 



(9coul ~ 2q cou iq a - 2Ag a )(g 2 oul 



2q CO uiq c — 2Ag c 



We observe that (for A = 0) this parameter £(0) is found 
to be negative and — £(0) ^> 1 for beam energy large com- 
pared to the binding energy and for perpendicular mo- 
mentum transfers q± 3> 2n a q^ (nonadiabatic case), where 



<7H = uj/v with Hu> = -E^ind + E rc i and where the relative 

H 2 a 2 

energy between c and n is E Te \ — with the relative 



momentum given by q 



It was already no- 



ticed in the numerical evaluation of the process that, due 
to — C(0) 3> 1 that the hypergeometric series does not con- 
verge and an analytic continuation had to be used. Here 
we use this fact to our advantage and make a linear trans- 
formation to get the argument of the hypergeometric func- 
tion close to 0. The transformation we are using leads to 
the argument of the hypergeometric function z = i_g(o) 
(Eq. 15.3.7 of Q). In this respect our approach differs 
from the one used in the bremsstrahlung case, where a 
transformation giving an argument close to one is used. 
Using only the lowest order term in the hypergeometric 
series one obtains after some algebra (up to an overall 
phase) 



T«47r A) /coui e~** 



qI - (q« + q c 



(8) 



+e 



m a q^ - (q« - q a ) 2 



Hereby, the relative phase is (f> — <To(t) c ) ~ cro(rj a ) — 0o(£) — 
f/21og|C(0)|. The 0-0(77) = argT(l + irj) are the usual 
Coulomb phase shifts, and £ = r) c — rj a . The correspon- 
dence to the Born result is clearly seen. One only has an 
additional prefactor and a relative phase e ±J * be- 

tween the two terms. The phase 4> obviously is 0(£). Since 
v c ~ v a the quantity £ is usually very small and so is <fi 
for the cases of 0,12). The prefactor is also well known 
in the semiclassical theory, where it accounts for the re- 
placement of the "Coulomb bended" trajectories with the 
straight line trajectories. Both corrections vanish in the 
limit £ — > and the result coincides with the usual Born 
approximation (even if rj a and rj c are not smalt). 

We have seen that the T-matrix in the case of large 
Coulomb parameters r\ a and r\ c corresponds to the Born 
result (small Coulomb parameter) in the sudden (or nona- 
diabatic) case q± 3> 2r] a q\\ . We note that the derivation of 
Eq. (||) only depends on the condition — £(0) 3> 1 (and not 
on the values of the r?'s). For n a ,rj c ^> 1 one can define a 
classical path for both a and c in the initial and final state 
and Eq. (||) can be related to the semiclassical approach 
(see the discussion below following Eq. (0)). We expect to 
find a connection between the semiclassical theory and the 
adiabatic case (q± < 2r} a q\\) and the fully quantal expres- 
sion for the T-matrix. For the adiabatic regime the well 
known exponential decrease with the adiabacity parame- 
ter is observed in the numerical calculations. In this case, 
the inequality — C(0) ^ 1 is not generally satisfied. We 
are presently working to see how the semiclassical limit 
can be obtained with analytical methods in this case also. 
Such a method would also be valid in both the adiabatic 
and nonadiabatic case as long as rj a , r\ c 3> 1 . 

A similar situation is encountered in the theory of 
bremsstrahlung and Coulomb excitation, see Section II E 
of (TJJ. There a fully quantal expression for the differen- 
tial cross-section for dipole Coulomb excitation is given in 
II E.62. It looks similar to the corresponding expression 
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for Coulomb breakup (see j| ) . The semiclassical " variant" 
of this formula is found in II E.57. It is noted there that 
it can be obtained from the quantal expression by letting 
at the same time r\ a and rj c go to infinity and perform a 
confluence in the hypergeometric functions. 



3 Scaling Properties 

In many experimental situations the Coulomb push q cou i 
is small. Having found that the full CWBA results agrees 
with the Born result for small scattering angles, we can 
expand Eq. (||) or Eq. (||) with <f> = £ = for small values 
of 5coui- We obtain 



coul " 



2D m^rric 2q • q c 

oul 

ml ( K 2 + (? 2)2 



(9) 



This result is in remarkable agreement with the usual I s * 
order treatment of electromagnetic excitation in the semi- 
classical approximation. 

In the semiclassical approach the scattering amplitude 
is given by the elastic scattering (Rutherford) amplitude 
times an excitation amplitude a(b), where the impact pa- 
rameter is related to the q± and r\, see above. The absolute 
square of a{b) gives the breakup probability P(b), in lowest 
order (LO). It is given by 



dPLo _ W 
dq 



3ttk (l+x 2 ) 4 



(10) 



where the variable x is related to the relative momentum 
between n and c by x = - and y is a strength parameter 
given by 

2ZZ c m n e 2 

y = ~t. r~ ■ ( n > 

nv a m a OK 

This formula shows very interesting scaling properties: 
Very many experiments, for neutron halo nuclei with dif- 
ferent binding energy, beam energy, scattering angles (or 
q„ and q c ) all lie on the same universal curvel (Correc- 
tions for finite values of £ e ff = tub/v = £(8) = 2n a q\\/q± 
should also be applied, according to E3].) It will be inter- 
esting to see in future calculations under what conditions 
(beam energy, . . . ) one finds deviations from this simple 
scaling behavior. E.g., postacceleration effects will lead to 
such scaling violations. 



4 Conclusion and Outlook 

The present model can be seen as a "theoretical labora- 
tory" , which allows to study analytically, as well as, nu- 
merically the relation between quantal and semiclassical 
theories, and the importance of postacceleration effects. 
We mention that from an experimental point of view, the 
postacceleration effects are not fully clarified, see, ej^ jll], 
|20| ("postacceleration") and on the other hand ("no 
postacceleration" ) . Finally, let us mention recent work on 
the electromagnetic dissociation of unstable neutron-rich 



oxygen isotopes |22|,|23|]. These authors deduce photoneu- 
tron cross-sections from their dissociation measurements. 
If the neutrons are emitted in a slow evaporation process 
in a later stage of the reaction, the question of postac- 
celeration is not there. On the other hand, for the light 
nuclei there is some direct neutron emission component 
and the present kind of theoretical analysis further proves 
the validity of the semiclassical approach used in p2| . 

Postacceleration effects are also of importance for the 
use of Coulomb dissociation for the study of radiative cap- 
ture reactions of astrophysical interest. We expect that 
our present investigations will shed light on questions of 
postacceleration and higher order effects in these cases 
also. 
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